An Almost L-Stable BDF-type Method for the Numerical Solution of Stiff ODEs Arising from the Method of Lines
نویسندگان
چکیده
A new BDF-type scheme is proposed for the numerical integration of the system of ordinary differential equations that arises in the Method of Lines solution of time-dependent partial differential equations. This system is usually stiff, so it is desirable for the numerical method to solve it to have good properties concerning stability. The method proposed in this article is almost L-stable and of algebraic order three. Numerical experiments illustrate the performance of the new method on different stiff systems of ODEs after discretizing in the space variable some PDE problems. © 2007Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 23: 1110–1121, 2007
منابع مشابه
On second derivative 3-stage Hermite--Birkhoff--Obrechkoff methods for stiff ODEs: A-stable up to order 10 with variable stepsize
Variable-step (VS) second derivative $k$-step $3$-stage Hermite--Birkhoff--Obrechkoff (HBO) methods of order $p=(k+3)$, denoted by HBO$(p)$ are constructed as a combination of linear $k$-step methods of order $(p-2)$ and a second derivative two-step diagonally implicit $3$-stage Hermite--Birkhoff method of order 5 (DIHB5) for solving stiff ordinary differential equations. The main reason for co...
متن کاملAn efficient method for the numerical solution of Helmholtz type general two point boundary value problems in ODEs
In this article, we propose and analyze a computational method for numerical solution of general two point boundary value problems. Method is tested on problems to ensure the computational eciency. We have compared numerical results with results obtained by other method in literature. We conclude that propose method is computationally ecient and eective.
متن کاملImplicit One-step L-stable Generalized Hybrid Methods for the Numerical Solution of First Order Initial Value problems
In this paper, we introduce the new class of implicit L-stable generalized hybrid methods for the numerical solution of first order initial value problems. We generalize the hybrid methods with utilize ynv directly in the right hand side of classical hybrid methods. The numerical experimentation showed that our method is considerably more efficient compared to well known methods used for the n...
متن کاملThe integration of stiff systems of ODEs using multistep methods
The second order Ordinary Differential Equation (ODE) system obtained after semidiscretizing the wavetype Partial Differential Equation (PDE) with the Finite Element Method (FEM), shows strong numerical stiffness. Although it can be integrated using Matlab ode-solvers, the function ode15s offered by Matlab for solving stiff ODE systems does not result very efficient as its resolution requires t...
متن کاملVariable Step Block Backward Differentiation Formula for Solving First Order Stiff ODEs
block method based on Backward Differentiation Formulae (BDF) of variable step size for solving first order stiff initial value problems (IVPs) for Ordinary Differential Equations (ODEs). In a 2-point Block Backward Differentiation Formula (BBDF), two solution values are produced simultaneously. Plots of their regions of absolute stability for the method are also presented. The efficiency of th...
متن کامل